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Many biological, ecological and economic systems are best described by weighted networks, as 
the nodes interact with each other with varying strength. However, most network models studied 
so far are binary, the link strength being either or 1. In this paper we introduce and investigate 
the scaling properties of a class of models which assign weights to the links as the network evolves. 
The combined numerical and analytical approach indicates that asymptotically the total weight 
distribution converges to the scaling behavior of the connectivity distribution, but this convergence 
is hampered by strong logarithmic corrections. 

PACS number: 84.35. +i, 05.40.-a, 02.50.Cw, 87.23. Ge 



With the increased availability of detailed topological 
data on complex networks describing communication, so- 
cial or biological systems emerges the need to develop 
tools to understand in general terms the origin and evo- 
lution of these complex evolving systems. The recogni- 
tion that real networks are fundamentally different from 
the random models that dominated the mathematical lit- 
erature in the past forty years lead to a surge of 
activity in addressing the statistical properties of these 
systems P-[To| . Despite rapid advances in uncovering the 
topology of complex networks, in one aspect most models 
are incomplete when compared with real systems: cur- 
rent network models assume that all links are equivalent. 
But in many fields it is well known that the interaction 
strengths can vary widely, such variations being essen- 
tial to the network's ability to carry on its basic func- 
tions. Sociologists have repeatedly argued about the im- 
portance of assigning strengths to social links, finding 
that the weak links people have outside their close circle 
of friends play a key role in keeping the social system 
together |11|. Recently, Newman has showed that as- 
signing weights to the links between scientists allows for 
a better characterization of the scientific collaboration 
web . Similarly, there is an ongoing discussion about 
the importance of weak links between species in guar- 
anteeing the stability of an ecosystem |13|] ; in addition 
viewing the economy as a complex network of companies, 
the monetary value of the various business transactions 
linking them is a key component of the systemic charac- 
terization. Finally, many transportation networks nature 
designed, ranging from cardiovascular to respiratory net- 
works, have well defined weights or flow rates assigned 
to the links, whose magnitude is intimately determined 
by the network's topology iQ. Despite the known im- 
portance of interaction strengths in various well known 
systems, there have been no attempts to model networks 
other than binary nets, whose links have weights or 1. 

In this paper we take a first step in the direction of 



a systematic study of networks with non-binary connec- 
tivities. We introduce and investigate two models that 
assign weights to new links as they are dynamically cre- 
ated, providing a prototype of a weighted evolving net- 
work. While we choose the simplest possible models, in 
which the weights are driven by the network connectivity 
only, numerical simulations indicate that the distribution 
of the total weight scales differently from the total con- 
nectivity. However, an analytical solution reveals that 
the different scaling behavior can be explained by strong 
logarithmic correction, and asymptotically the investi- 
gated weighted networks belong to the same universality 
class as their unweighted counterparts. 

Weighted Scale-Free (WSF) Model : Starting from a 
small number (toq) of vertices, at each time step we add 
a new node which links to m existing nodes in the sys- 
tem. The probability that a new node j will connect to 
a existing node i is 



ki 



(1) 



where ki is the total number of links that the node i has. 
In assigning a weight to the newly established link j ^ i, 
we assume that the weight Wji{= Wij) is proportional to 
ki, i.e., more connected (and therefore more "powerful") 
nodes gain more weight. Also, one can assume that all 
new nodes have fairly uniform total 'resources' for linking 
to other nodes in the system, we therefore require that 
each new node has a fixed total weight, i.e. we normalize 
Wij such that the sum of the weights for the m new links 



= 1, where {«'} represents a sum over the 



existing nodes to which the new node j is connected. As 
a result of the two assumptions, each link i <-> j of the 
newly added node j is assigned a weight as 



E 



(2) 



Weighted Exponential (WE) Model : The model is in- 
spired by model A discussed in Refs. jl^Jl^, and is de- 
fined as follows : at every time step we add a new node 
with m(< Too) links, connected with equal probability to 
the nodes present in the system. The weights of the links 
are assigned again by using (||). 

The difference between the WSF and WE models 
comes in preferential attachment, which is known to fun- 
damentally alter the topology : The WSF 
model generates a scale-free network whose connectiv- 
ity distribution follows P{k) ^ k"^, while the network 
generated by the WE model is exponential with the con- 
nectivity distribution following P{k) — ■^e~^/™. Since 
the weights of the links are driven by the connectivity, 
this difference is expected to lead to significant changes 
in the distribution of the link strengths as well. 

We start by investigating the weight distribution of the 
two models. As Fig 1 (a) and (b) shows, both the WE 
and the WSF models lead to a peaked and skewed weight 
distribution, whose tails decay exponentially (or faster) 
for large wtj. The boundedness of P{wij) is due to the 
normalization condition, which does not allow individual 
weights to be larger than 1. Most important, however, 
we find that the distribution is stationary, i.e. P(wij) is 
independent of time (and system size). 

While the individual weights assigned to links, Wij, 
are bounded, we get a very different picture when we 
study the total weight associated with a selected node. 
In binary networks the node's importance is character- 
ized by the total number of links it has, ki. Similarly, in 
a weighted network the importance of a node i can be 
measured by its total weight, obtained by summing the 
weights of the links that connect to it, Wi = ""^ij- 

Due to the normalization condition (^) a new node has 
Wi — 1, but Wi increases in time every time when a sub- 
sequently added nodes link to i. Since in both models 
the weights are determined by the network connectivity, 
we expect that P{w) closely follows P{k). In contrast, 
the numerical results summarized in Fig. 2 indicate strik- 
ing differences between P{k) and P{w). As Fig 2a shows, 
while for the WE model P{k) decays exponentially, P{w) 
systematically deviates from a simple exponential behav- 
ior. This difference is even more evident in the network 
dynamics: while both ki{t) and Wi(t) appear to increase 
logarithmically in time, they can be fitted with a different 
slope on a log- linear plot (Fig. 2b). Similar systematic 
discrepancies are observed for the WSF model as well: as 
Fig. 2c indicates, while P{w) can be fitted by a power 
law, P{w) ~ w^'^ , it appears that the exponent a is dif- 
ferent from 7 = 3. Furthermore, we find that a depends 
strongly on to (Fig. 2c). Again, this difference is re- 
fiected in the dynamical behavior of ki{t) and Wi{t): as 
Fig. 2e indicates, Wjjt) t^ with /? > 1/2, in contrast 
with ki{t) ^ t^^^ predicted by the binary scale- free 

model. 

To understand the different behaviors of Wi and ki un- 
covered by the numerical simulations, we resort to ana- 



lytical method in determining the averaged behavior of 
Wi{t) for the discussed model. To simplify the discussion 
in the following we assume m = 2, however, the calcu- 
lations can be generalized for arbitrary to. The total 
weight of node i at time t can be written as 



it) = l + J2^v A it') {w,, it')) dt', (3) 



where Pi{t) is the probability that node i is selected to 
be connected to a new node j at time t and is the 
time at which the node i has been added to the system. 
{wij ) is the average weight of link i ^ j once the link is 
established. When a new node j and the list of to nodes 
{i'} to which it connects are selected, the weights of the 
links, Wji' are assigned according to (^). These weights 
depend on the number of links the selected nodes have, 
i.e. {fci/}. If we assume that node j is connected to nodes 
i and I (to = 2), we have 



(wy(t)) = / w,,{l)V{ki) dki 



(4) 



where Wji(l) is the weight between the j and i nodes, 
Viki) is the probability distribution of ki, the total link 
number of node I. Substituting (0) into (0), we obtain 



Wi{t) 



m')wj,{l)V{ki) dki dt'. (5) 



by 



According to (0) for to = 2, the weight Wji{l) is given 



Wji{l) 



thus Eq.(||) becomes 



W^{t) = 1 + 



k 



(6) 



■V{ki) dki dt'. (7) 



Eq. (|^) represents a general expression for calculating 
Wi{t) for TO = 2. To apply it to the WE and WSF models, 
we need to calculate explicitly P{t) and V{ki). 

WE model: In the WE model the nodes to which a 
new node connects to are selected uniformly among all 
existing nodes, thus the probability that node i will be 
picked is independent of this node's connectivity and is 
given by 



my- 



t + niQ 



(8) 



Similarly, the connectivity distribution and the dynami- 
cal behavior of a single node are given by 



e 

TO 



fc/m^ (9) 



ki{t)= TO [ln(TOo + i - 1) - ln(TOo + t" - 1) + l] 
= 771 [ln(ai + 6) + 1] , 



2 



where a 



mo+t"-l • 



b = u"'" , and the normahzation 



Wi 



condition is 1 — V{k)dk. 

Substituting (|) into (|), we obtain 



w.,{t) = 1 + e 



hit') 



t« J,n t' + mo ki{t') + ki 



6"'='/'" dk, dt'. 



After performing the integration and inserting ki (t) from 
for large t we obtain 

Wi{t) ~mln(at + 6)-mln(ln(ai + &) + 2) + C, (10) 

where C is an integration constant independent of t. 
Therefore the relation between Wi{t) and ki{t) for large 
t follows 



Wi (t) ~ ki (t) — ni In In t + C. 



(11) 



The prediction ( |ll| ) is fully supported by numerical sim- 
ulations: in Fig. 3a we plot the difference Wi{t) — ki{t) 
as function of Inln(t), showing that the difference indeed 
follows a double logarithmic law. This result is very in- 
teresting since it indicates that the different slopes ob- 
served in Fig. 2b for ki{t) and Wi{t) do not represent 
distinct power law scaling behaviors, but are the result 
of logarithmic corrections. 

WSF model: In the scale-free model the probability 
distributions and ki (t) are given by 

~ _fc^ hit) hit) 
Fi{t)= m — ^ = TO- — 



2mt 2t 



e; h 

Vik)^mk-^ (cxfc-P(fc)) 



(12) 



kdt)-- 



,Vt. 



Substituting ^ into (0), and performing the integrals 
we obtain 



TO / , TO t \ TO 



W^it) ~ hit) - — (In 



+ _lnmln^ + C", (13) 



indicating that despite a different scaling behavior sug- 
gested by the numerical simulations (Fig. 2e), we are 
dealing with strong logarithmic corrections and asymp- 
totically we have (3' — (3. Again, the analytical prediction 
( p^ ) is confirmed by more detailed numerical simulations 
shown in Fig. 3b. 

Our ability to calculate analytically Wij for the dis- 
cussed models is based on the fact that the weights are 
driven by the connectivity distribution. To address the 
generality of our results we investigated several exten- 
sions of these two models, that we discuss in the follow- 
ing. 

Weight driven weight- In general one could expect that 
in some systems the quantity determining the weight is 
not the connectivity, but are the weights themselves. To 
investigate this possibility we replaced (S) with 



(14) 



i.e. the weight of the newly added links are determined 
by the total weight of the nodes. While we cannot solve 
this model analytically, the numerical results are simi- 
lar to those observed for the WE and WSF models: an 
apparently different scaling behavior for k and w can be 
attributed to slow corrections to scaling. 

Weight driven connectivity- In some systems the topol- 
ogy could be driven by the total weights, and not by the 
connectivity. Thus we assume that the probability (0) 
that a new node is connected to a node j is 



Wi 



Ej ^0 



(15) 



where Wi is the weight of node i. The weights are then 
assigned following 1^. We find that the scaling of this 
network is identical to that of the scale-free model, and 
the evolution of the weights also follows the paradigm 
established for the WSF model. 

Discussion- Weighted links are a common feature of 
real networks, thus addressing their scaling behavior is 
of primarily importance if we are to understand com- 
plex networks in general. Here we take a first step in 
this direction by investigating the scaling properties of 
several simple models that incorporate mechanisms to 
assign weights to the links. The weight distribution in 
the discussed models can be determined analytically and 
by using numerical simulations, allowing for a thorough 
analysis of the model's scaling properties. A compari- 
son between the analytical and numerical results brings 
deeper understanding of the observed scaling behavior 
for the distribution functions. First, extensive simula- 
tions of networks whose size is comparable to the real 
networks that are currently available indicate the emer- 
gence of new scaling exponents for the behavior of the 
total weights. However, the analytical solutions reveal 
that the results are affected by strong logarithmic cor- 
rections, and asymptotically the scaling behaviors of the 
weighted and unweighted models are identical. This re- 
sult raises important questions regarding our ability to 
uncover the correct scaling behavior of real weighted net- 
works, should such data become available in the near 
future: the real exponents could be easily shadowed by 
corrections to scaling similar to that encountered in the 
investigated models here. 

An important feature of the studied models is the fact 
that the distribution of the individual weights, w^ , is 
bounded and stationary. This implies that the scaling 
behavior of the total weight distribution can not be ex- 
plained by the distribution of individual weights. Indeed 
it is the correlation between the total weight and the 
linkage probability of the new node that drives the dis- 
tribution towards power law behavior. In general, one 
can imagine systems and models where the link distri- 
bution, Piwij), could also have nontrivial distribution. 
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such as a power law. Investigating such models could be 
of major future interest. 

The results presented in this paper represent only the 
starting point towards understanding weighted networks. 
In some real systems, diverse dynamical rules can govern 
the assignment of weights to links, which could result 
in statistical properties of the network that are differ- 
ent from that discussed here. In particular, we assumed 
that once a weight has been assigned to a link, it stays 
unchanged, which is often not the case in more realis- 
tic networks: weights can evolve dynamically just as the 
network topology does. For example, acquaintance can 
turn into friendship by strengthening a previously weak 
link. Determining the generic behavior of such complex 
evolving systems is a real challenge for future research. 
Despite these limitations, the investigated models give a 
glimpse into the complex behavior we are facing as we 
attempt to make network modeling more realistic by in- 
corporating weights. Finally, the ultimate understanding 
of weighted networks will be determined by the available 
data on real systems. While currently such data is rare, 
we believe that the increasing interest in network model- 
ing and creative data collection methods will soon lead to 
the development of such data sets, offering further guid- 
ance for modeling these complex systems. 

Research at Notre Dame supported by NSF, PHY- 
9988674 and CAREER DMR97-01998. 
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FIG. 1. The distribution P{wij) of the individual link 
weights, for the (a) WE and the (b) WSF models, defined 
in the text (m — 2). The symbols correspond to different 
system sizes (or time), i.e. = lO^(O), W^i^), W^iO) and 
10^ (A). The insets shows the same data on a log —linear plot, 
indicating that the tail decays faster than exponential. 
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FIG. 2. (a) Distribution P{w) of the total connectivity w 
assigned to individual nodes for the WE model. The symbols 
correspond to different values of m, i.e. m = 2 (O); 3 (^)i 
4 (0) and 5 (A). The inset shows the connectivity distribu- 
tion, P{k), for the same parameters as in the main panel, (b) 
Time dependence of ki{t) (Q) and Wi{t) (□) for a randomly 
selected node i for the WE model {i = 5000). (c)P(fc) (Q) 
and P{w) (□) distributions for the WSF model for m = 5. 
The inset shows the same data for m = 2. (d) ki{t) (0)i 
w^{t) (□) VS. t for the WSF model {i = 10000). 



4 





WSF 


(b) 





FIG. 3. The diflference {wi{t) - hit)) for the (a) WE and 
the (b) WSF models. The continuous lines in each case repre- 
sent the analytic solution (11) and (13), respectively. Wc lim- 
ited the simulations to nodes appears at large t^it^l = 10*) to 
capture the asymptotic limit, that is predicted by our predic- 
tions (11) and (13). We find that for smaller f° the crossover 
time for the convergence to the analytic solution is numeri- 
cally prohibited. 
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